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Abstract 

The Sznajd model is a sociophysics model that mimics the propagation of opin- 
ions in a closed society, where the interactions favour groups of agreeing people. 
It is based in the Ising and Potts ferromagnetic models and although the original 
model used only linear chains, it has since been adapted to general networks. This 
model has a very rich transient, that has been used to model several aspects of 
elections, but its stationary states are always consensus states. In order to model 
more complex behaviours we have, in a recent work, introduced the idea of biases 
and prejudices to the Sznajd model, by generalizing the bounded confidence rule 
that is common to many continuous opinion models. In that work we have found 
that the mean-field version of this model (corresponding to a complete network) 
allows for stationary states where non-interacting opinions survive, but never for 
the coexistence of interacting opinions. In the present work, we provide networks 
that allow for the coexistence of interacting opinions. Moreover, we show that 
the model does not become inactive, that is, the opinions keep changing, even in 
the stationary regime. We also provide results that give some insights on how this 
behaviour approaches the mean-field behaviour, as the networks are changed. 

1 Introduction 

In the last decade, the study of complex networks and of opinion propagation models 
has become increasingly common. Opinion propagation models seek to explain prop- 
erties of elections, the spread of rumours and the formation of factions and consensuses 
in communities. These models, dubbed sociophysics models, apply approaches typical 
of statistical mechanics, using Monte Carlo simulations and microscopical rules. 

The Sznajd model was defined by Sznajd and Sznajd- Weron under the name "United 
We Stand. Divided We Fall" as a model where the society is represented by a linear 
chain, and people can have one of two opposite opinions |4 |. This model was quickly 
adapted to an arbitrary number of opinions and a general network J2). The main dif- 
ferences between the Sznajd model and other opinion propagation models are that the 
interactions favour bigger groups of agreeing people and the agents can be seen to 
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be influencing their environment, instead of the opposite (which happens in the voter 
model, for instance). 

The transient of this model shows some interesting properties, specially when the 
society is modelled by a complex network, and has been used to explain some scaling 
properties of both proportional and majority elections [2j . On the other hand, the sta- 
tionary states of the model are always consensus states, implying that everyone has the 
same opinion. 

In a recent work 0, we have generalized the rules of the model, in a way that 
can be interpreted as the introduction of prejudices and biases. We found out that 
the mean-field behaviour of this generalized model allows for the coexistence of non- 
interacting opinions in the stationary state, and that societies modelled by Barabasi- 
Albert networks [ 1 1 follow essentialy the mean-field behaviour. 

In the present work we show that some Watts-Strogatz networks |5), based on 
square networks, allow for stationary regimes, where not only interacting opinions 
coexist, but the quantity of nodes holding a given opinion oscillates with time. 

2 Model Description 

We now introduce the Sznajd model as defined in @, the bounded confidence rule IT) , 
general confidence rules and a special case that can be interpreted as 2 identical 
parties, that have each 2 factions with different convincing powers. 

2.1 The Sznajd Model 

In the Sznajd model the society is represented by a network, where each node in the 
network is a person, each edge is a social connection (friendship, acquaintance, mar- 
riage, etc.) and each node i is associated to an integer variable er,;, between 1 and M, 
representing the opinion of the corresponding person. 

At each time step, we choose a node i at random and (also at random) a neighbour 
j of i. If i and j disagree (<Tj ^ cry), nothing happens. Otherwise, we choose at 
random a neighbour k of j and the node k becomes convinced of the opinion of the 
pair ((Tfc — > Gj) with some probability p. In the usual model, p is always 1. Figure Q] 
shows a scheme of the update rule in action. 




Figure 1: The Sznajd update scheme. First a node i is chosen at random. Then a 
neighbour of i, j, is chosen. Finally, if i and j have the same opinion, they attempt to 
convince a neighbour of j, k, chosen at random. 



2 



2.2 Bounded Confidence and Confidence Rules 



The bounded confidence rule was introduced in (3] as a way to tell how close 2 con- 
tinuous opinions are to each other and to model extremism. This rule says that opinion 
changes are not abrupt, which is accomplished by imposing that 2 agents with opinions 
cr and a' only interact if \cr — a'\ < e. This rule can be adapted to discrete opinion 
models, like the Sznajd model, by taking e = 1 and changing the number of opinions 
instead of the interaction threshold. In the Sznajd model this means that a can only be 
convinced to an opinion a' iff a' = a ± 1. 

In 0, we have generalized the bounded confidence rule to what we called gen- 
eral confidence rules, by imposing that the probability that an opinion a be convinced 
to an opinion a' is given by p a ^a' (a parameter of the model). These rules can be 
schematized by interpreting these parameters as the elements of the adjacency matrix 
of a weighted directed graph, and associating this graph to the set of parameters (this 
graph will also be refered to as the confidence rule). 

It is noteworthy that many of the previous modifications of the Sznajd Model (as 
well as the usual version, introduced in section 12. It can be obtained by appropriate 
choices of the parameters: 

• The usual model is the case p a >^a = 1 for all a ^ a' . 

• If an opinion a has p a '^a = and p a ^a' = P 7^ V a', a behaves analogously 
to an undecided state. 

• Usual bounded confidence is the case p a >^ a = lif<x' = er±l and p a '^ a = 
otherwise. 

We also introduced a phase space representation, in which the phase space variables 
are the proportions of sites in the network with a given opinion. We will denote these 
variables r) a , that is, r\ a is the proportion of sites in the network with opinion cr. These 
variables have the following constraints: 

y^?7ff = 1 and J fr > 0- 

cr 

which imply that the phase space must be a simplex (the phase space with 2 opinions 
is a line segment, with 3 opinions it is a triangle, a tetraedron for 4 opinions and so on). 

In the next section we investigate a family of confidence rules that can be inter- 
preted as a dispute among 2 identical parties, that have 2 different factions each that 
behave differently. This family is characterized by 3 parameters, p, q and r, and is 
given by: 

Pl->2 = P2^3 = P3^4 = Pi^l = P 
P2^1 = PZ^I = Pi^3 = Pl^4 = q 
Pl^3 = P3^1 = P2^A = PA^2 = V. 

Figure [2] illustrates this rule, and the 2 parties {1,3} and {2,4}, as well as the 
symmetry between them. We will denote this rule R(p,q,r). The mean-field re- 
sults say that if r = we have 2 attractors: {(771, 772, f?3, V^lvi + V3 = 1} ar, d 
{('7i, ^2, ?73,?74) I ??2+?74 = l};ifr ^ Owe have 4 consensus attractors ({(771, 772, 773, 774) 1 77^ = 
1} for a = 1, 2, 3, 4); and if p = q = we have a degenerate case, with 2 independent 
models following the original confidence rule (we will ignore this case). 
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Figure 2: The 2 parties rule. Note that this graph has rotation symmetry, hence the 
2 parties (green and red) are identical. The mean field results state that one of the 2 
parties will always dominate the whole network and unless r — 0, only one of the 
factions survive. 



3 Simulations 

We use for these simulations a network based in the Watts-Strogatz model J6), but in- 
stead of rewiring a ring like network, we start from a square network with periodic 
boundary conditions and rewire s.E edges, chosen at random (s is the rewiring param- 
eter and E is the number of edges). All the networks used have, as a starting point, a 
square network with side L = 316. 

For the initial conditions of the simulations, instead of drawing the opinions uni- 
formly, we follow the approach introduced in [5] . Firstly, we draw a point in the tetrae- 
dron {(771,772, m, f74)|77i+772 +773 + 774 = 1 and 771, 772, V3, m > 0} uniformly and then 
we draw the opinions 1, 2, 3 and 4 with probabilities 771, 772, 773 and 774 respectively. 

We found that while some of the simulations end up in an inactive stationary 
regime, an oscillating (and therefore, active) stationary regime is also possible. Fig- 
ure [3] shows the time series for the supporters of one of the 2 parties, using a fixed 
confidence rule, where this behaviour can be seen. 

The confidence rule depends on 3 parameters, however the stationary qualitative 
properties can be studied in only 2 dimensions, because the only difference between 
R(p, q, r) and R(Xp, Xq, Xr) is the time scale in which things happen. Also R(p, q, r) 
and R(q,p, r) are actually the same rule, as the corresponding graphs are isomorphic 
(you can go from one rule to the other by swapping opinions 2 and 4, for example). Fig- 
ure [4] shows the percentage of simulations that enter the oscillating stationary regime 
for some values of ^ and a (the rewiring parameter s was kept fixed at s = 3.10 -3 ). 
It shows that these oscillations vanish when p ~ q, meaning that some degree of as- 
symetry is necessary in the confidence rule for this regime to appear. 

We now keep the confidence rule fixed (we use p = 1, q = 0.1 and r = 0) and 
change the rewiring parameter. Also, as the rewiring parameter grows, the average 
path lenght decreases and new parts of the network become connected. As we in- 
crease rewiring, the amplitude of the oscillations (which we defined as the average of 
the moving standard deviations over the oscillating simulations) grows, less simula- 
tions reach the oscillating stationary regime and eventually the regime ceases. Figure 
[5] shows the percentage of simulations that reach an oscillating stationary regime and 
their amplitude as a function of the average path lenght I. 

Figure [3] shows the typical behaviour of the oscillations in a case with higher aver- 
age path lenght (i ~ 42.6) and figure |7]for a case closer to the transition to the non- 
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Figure 3: Time series for rji + 773 in the case p = 1, q = 0.1 and r = 0. The rewiring 
parameter is s = 10~ 3 and the network's average path lenght is approximately 42.6. 
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Figure 4: The percentage of simulations that end up in the oscillating stationary regime 
as a function of £ and 2 . The network is fixed, with rewiring parameter s = 3.10 _3 and 
average path lenght approximately 29.8. The graph was obtained using cubic splines. 
The marks in the graph represent the actual data points. 
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Figure 5: Amplitude of oscillations and percentage of simulations ending in the os- 
cillating stationary regime. Note that there where no oscillating simulations found for 
I ~ 14.8 and I ~ 11.0, corresponding to rewiring parameters s = 0.03 and s = 0.1 
respectively. The gray region is where the transition from existence to abscence of the 
regime is observed. 
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oscillatory regime (I ~ 20.3). These results can be understood from a dynamical sys- 
tems perspective, using { (7/1 , r? 2 , m , ij 4 ) I m + m + m + m = 1 and m , m , m , m > o } 
as the phase space. We will use projections of this 3 dimensional phase space (a tetrae- 
dron) into 2 dimensions, in order to get a clearer picture of the trajectories, depicted in 
those time series. The phase space we are considering is a simplex, so each point in it 
is a convex linear combination of its vertices. Let Pi be the vertices of the phase space, 
let Q be a projection of this phase space (a linear mapping) into some vector space and 
let Qi = QPi. We can choose to embed our phase space in a way such that the point 
corresponding to (771 , 772, 773, 774) is 

a 

So if we project the phase space, we get 

a a 

So we can either choose the linear mapping Q or the coordinates of the projected 
vertices, Q a . To build our projection of the 3 dimensional phase space in R 2 we choose 
Qi = (0,1), Q2 = (1,0), Q3 = (0, — l)andQ4 = (—1,0), meaning that the projected 
phase space is a square and the time series are obtained plotting r\\ — 773 against 772 — 774. 
Geometricaly, this corresponds to the projection depicted in figure|6] 




Figure 6: Scheme depicting the projection that will be used to visualize the trajectories 

The oscillations in figures [3] and [7] are due to an attractor resembling a limit cycle, 
whose size increases as the rewiring parameter is increased (figures [8] [9] and [TOt . but 
whose attraction basin becomes smaller (these 2 conditions can happen at the same 
time, if the basin is becoming thinner). 

Finaly we present the result of these projections for a Barabasi- Albert network 
ID, as a comparison (figure [TTV We can see trajectories connecting the attractors 
(the vertices of the square) which are reminiscent to the attractor responsible for the 
oscillating stationary regime. 
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Figure 7: Time series for r/i + 7/3 in the case p = 1, q = 0.1 and r = 0. The rewiring 
parameter is s = 10~ 2 and the network's average path lenght is approximately 20.3. 




Figure 8: Projection of the trajectories in phase space for the network with s = 10~ 3 
and I ~ 42.6. The projections of the mean field attractors are the line connecting the 
red points (coexistence of opinions 1 and 3) and the one connecting the green points 
(coexistence of opinions 2 and 4). 
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Projection of the trajectories 




Figure 9: Projection of the trajectories in phase space for the network with s = 3.10 -3 
and £ ~ 29.8. The projections of the mean field attractors are the line connecting the 
red points (coexistence of opinions 1 and 3) and the one connecting the green points 
(coexistence of opinions 2 and 4). 
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Projection of the trajectories 




Figure 10: Projection of the trajectories in phase space for the network with s = 10 
and I ~ 20.3. The projections of the mean field attractors are the line connecting the 
red points (coexistence of opinions 1 and 3) and the one connecting the green points 
(coexistence of opinions 2 and 4). 
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Projection of the trajectories 




Figure 11: Projection of the trajectories in phase space for a Barabasi- Albert network 
with minimal degree 5, network size 10 5 and I ~ 4.3. We use for this projection 
p = 0.125 and q = r = 1. The mean field attractors are the vertices of the square. 
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4 Conclusions 



We have shown that the Sznajd model, generalized to include prejudices and biases, 
can have active stationary states, with the proportion of sites holding a given opinion 
oscillating over time. This oscillatory regime seems to depend crucially on assymetries 
of the confidence rule used and on the average path lenght not being too small. The 
model simulated in the square lattice displays these oscillations and they remain, even 
if the edges are rewired (provided they are not rewired too much). The oscillations 
seem to be due to an attractor, resembling a limit cycle. The radius of the cycle in- 
creases as the rewiring parameter grows and the average path lenght decreases, but the 
volume of the basin of attraction becomes smaller, probably indicating that it is thin- 
ning. This seems to indicate that the transition between the existence and the absence 
of an oscillating stationary regime is due to either the limit cycle colliding with the 
mean-field attractors, or the basin of attraction collapsing (together with the attractor) 
to the frontier between the basins of the mean-field attractors (these hypothesis seem 
to be supported by figure [10] and to some extent by figure [TT). 
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